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Abstract: We investigate the non-Pauli-Fierz (nPF) theory, a linearized massive gravity 
with a generic graviton mass term, which has been ignored due to a ghost in its spectrum 
and the resultant loss of unitarity. We first show that it is possible to use the Lee- Wick 
mechanism, a unitarization through the decay of a ghost, in order to handle the sixth mode 
ghost of nPF, and then check for the quantum consistency. Once proven to be consistent, 
nPF could become a viable candidate for a large distance modification of gravity, because it 
naturally solves the intrinsic problems that most dark energy /modified gravity models suffer 
from: It smoothly converges to general relativity at short distances, and the small graviton 
mass necessary to modify gravity at large scales can be stable under the radiative corrections 
from the minimal gravity-to-matter coupling. 
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1 Introduction 

Although the Pauli-Fierz massive gravityfl] is special due to the absence of a ghost, a particle 
or field with a wrong-sign kinetic term, it is hard to keep it that way. While whether the 
Boulware-Deser ghost [2] can be eliminated at all nonlinear orders is still in debate [3], we 
saw[4] that quantum loop corrections revived the ghost for a large class of interactions. Then, 
instead of trying to obtain a ghost-free theory, how about accepting the presence of a ghost 
and finding a way to live with it? 

Allowing a ghost to appear in the spectrum, we get the non-Pauli-Fierz (nPF) massive 
gravity. Its action in a flat 4d background is 



with h = rj^h^v and rj^ u = diag(— 1, 1, 1, 1). The corresponding tree level propagator is 




(1.1) 




(1.2) 
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whose first line contains a massive spin-2 graviton, whereas the second line shows a spin-0 
ghost. Here Ij is a set of tensor bases with 4 indicies: 



Ii = f/> Vx P > h = VfiXVup + VwVu\ , I3 = r]p\k„k p + r] pp k u k x + (/i -B- z/) , 

I4 = Vpykxkp + k p k u rj\p , I5 = k p k v k\kp . (1-3) 

Unlike the usual spin-projection operator of, e.g., [5], Ij's are not orthogonal, but still complete 
because any symmetric second rank tensor can be expanded in terms of them. When coupled 
to a conserved source, the k S> m limit of (1.2) gives the massless graviton propagator, 
4j( — y + if)> °f general relativity (GR) in de Donger gauge, which signifies that nPF 
will agree with GR at short distances. To check this claim, let us repeat van Dam and 
Veltman's analysis, [6], for nPF: Consider a gravitational interaction between two conserved, 
non-relativistic sources, T and T 1 . With <7nPF the coupling of an nPF graviton to T and T', 
the tree level contribution to this process is 



1 2 T 00 T' 00 1 
ffnPF 73. ^ 17) \~U> ) " 1-4 



2 ^ fc2 . fc2 -VA;2 
Since the same process mediated by a massless graviton is 

I ^OOjvOO 

< L5 > 

where g is the gravitational coupling of GR, we see that <?nPF = 9, z.e., nPF matches GR at 
distances much larger than the inverse graviton mass. 1 Thus, if we manage to recover unitarity 
that is apparently violated by the ghost, nPF might become a viable candidate for a large 
distance modification of gravity: GR is restored at short scales, while the gravitational force 
gets weakened by Yukawa-suppression at distances 3> m _1 , so that long wavelength sources 
such as the cosmological constant gravitate less, providing a solution to the cosmological 
constant problem in the sense of [7]. 

Then, the spectrum of nPF suggests to pick up, among various proposals[8] to deal 
with ghosts, the Lee- Wick (LW) mechanism[9] as a resolution to its ghost problem. The LW 
mechanism unitarizes a theory with a ghost by the decay of the ghost into ordinary fields. 
Mainly as a means of dealing with ghosts from higher derivative theories, it has been applied 
to constructing a renormalizable quantum gravity[10] as well as a hierarchy problem-free 
standard modelfll]. To elucidate, let us highlight a scalar toy model worked out in [11]: 

C = -\dpad»a - \rn\a 2 + \d^p + ^m 2 pP 2 - |(a - pf , (1.6) 



1 In comparison, Pauli-Fierz massive gravity gives gp F — |g 2 , which is the source of the notorious van-Dam- 
Veltman-Zakharov discontinuity. 
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where o is an ordinary scalar and p is a ghost. Due to the cubic interaction, the p propagator 
gets renormalized by the self-energy, — zS, 



P M = -5 5 — ^T^T • (1-7) 

p p 2 + mj + T,(p 2 ) v ; 

When the p — > oo decay channel is open, i.e., m p > 2m a , a contribution to S from the 1-loop 
self-energy with two poo- vertices has an imaginary part, £ = im p T, 



Now let us consider, for example, the scattering amplitude, Ai, of oo — > oo scattering. There 
is a possibility of unitarity violation when the scattering is mediated by a virtual p: 

iM = {-ig) ■ P p ■ {-ig) . (1.9) 

But the imaginary part of Ai for this process 

— (p 1 TTi r 

ImM = , (1.10) 

{p A + m,pY + rripl z 

is still positive because the unconventional sign of T is compensated by the unconventional 
sign of the residue of P p , and therefore unitarity is maintained. A complete analysis of 
causality, stability, etc., as well as unitarity, of LW theory needs more subtle and elaborate 
considerations, and has been studied in the literature[12]-[13]. 

It seems possible that the same mechanism would work for nPF when 2(i-a) m2 ^ 4m 2 . 
In the next section, we will see our expectation being fulfilled. Then in §3, we check a basic 
quantum consistency of nPF. We discuss on the possible issues and future directions in §4. 
The details of the loop computation can be found in the Appendix. 

2 LW mechanism for nPF 

Let us first see if or how the LW mechanism works for nPF. We start with (1.1) and decompose 
hav in the usual way: 

Ku=fuu + — ~ + -^ b77/^02, (2.1) 

m m z 

with d^fftv = 0, / = Tf v fp U = and d^V 1 = 0. Then the onshell constraint 



.„<5SnPF 

u = c 

gives 



= ^7^f = d»h IJlv -ad IM h, (2.2) 



— - + o + d^<j)2 = ad^ — T + # 2 , (2.3) 
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or 

r,-o, ^ = ip± 02 . (2.4) 

m z 1 — a 

With (2.1) and (2.4), (1.1) becomes 

from which we identify the physical degrees of freedom(DOFs) to be a spin-2, /™, of mass 
m and a ghost scalar, fa, of mass m a = y ^i-a) m ' Then we introduce a generic quantum 
interaction 



1 d^AfCi^ 1 ^^ + C 2 V^ + Cs^ 3 ) , (2-6) 



to obtain 



s = J d 4 x{ - \r{-d 2 + m 2 )i 2 /^ + 3^ 2 ( - a 2 + ^j^y™ 2 )^ 

+A(3Ci + C 2 )/^r>2 + •••}, (2.7) 

where • • • is the interactions of the form of / 3 and (fy. With 

Pi I 2 „ — i 
£ — p — 

/ 9 I 9 o — ' 02 9 i 9 — ' 

pZ _|_ m z 2 77i^ 

V //0a = -iA(3Ci + C 2 )y = >-, 
we get, for A; 2 pa — m 2 , 

d 4 p r .., OA j. ,-,2l 2 i h I 2 « I2 



(2vr) 4 1 vs s ;J 2 p 2 + m 2 2 2 (p + k) 2 + m 2 2 

5a 2 (3 Cl + C2) 2 ( ^ + _Ji -_ + ...) , (2. 



m 2 . 



where we use dimensional regularization and • • • is a finite real contribution. In the narrow 
width approximation, (2.8) defines the decay width of fa into two /'s: 

—i —i 

P* 2 Pfc.rcnorm = p 2 + m 2 + £ = p 2 + m 2 + -^f + . . . > ( 2 ' 9 ) 

with 



r = _<^26±C01/ 1 _±^. (2 . 10) 

Therefore, nPF theory satisfies the conditions for the LW mechanism to work when m a > 2m, 
i.e., I < a < 1, and the imaginary part of the matrix element of the // — > ff scattering 
mediated by fa, 

i{-i\(3(i + C2)} 2 TT • ^2,rcnorm ■ Tt] « To~, 9TT~ 9~F9 ' i 2 ' 11 ) 



InxM = Im 

is positive, as required by unitarity. 



(p 2 + m 2 ) 2 + m 2 T 2 ' 
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3 Renormalizability 



Since quantum effects play an essential role in rendering nPF theory sensible, we have to check 
the quantum consistency of the theory. One of the items in the checklist is renormalizability. 
Divergences of the loop corrections to SnPF would have pieces of 0(k i ) and higher orders as 
well as O(k ) and 0(k 2 ) terms. Adopting the idea[14] of effective field theory (EFT), we 
may assume those divergences of higher orders in k can be absorbed into EFT terms. But 
regardless of the EFT treatment, O(k ) and 0(k 2 ) divergences must be taken care of at tree 
level. With the divergences of the 1PI 2-point function written as 

n div = -{&l, m Il + 2 , m I 2 + £ 2 (&l,fcll + b 2>k I 2 ) + 63I3 + 64I4 + o{k 4 )} , (3.1) 

the very first requirement for any desired interaction is that its 0(k 2 ) part conform to the 
tree level form, k 2 {l\ — ^) + — I4, i-e., 

h, k + 2b 2 ,k = 0, h >k + h = 0, 6 2>fc + 63 = . (3.2) 

When there is a 1PI 1-point function(tadpole), 

B d iv = -hv^ , (3.3) 



it should be interpreted as the 0(h) piece of Ay — det^j, + jffj-) = A(l + ^g^- + ^ 



4M <2 — h • • •)■ Then the conformity at O(k ) means 

b h + h- (b 1:m h + &2,ml 2 ) • h (3.4) 

should be able to be written as 



A 



for some A and X. That is, the second requirement is 

b 1>m + 2ab 2 , m = (1 - 2a )^^ ■ ( 3 - 6 ) 
3.1 Polynomial interaction 

At 1-loop, the 0(k 2 ) part of (3.1) for the cubic polynomial interaction, (2.6), is calculated in 
A.l, from which we can obtain the explicit expression for (3.2): 



= 9(37 - 74a + 52a 2 + 84a 3 )( 2 



1 

+4C 2 {(193 - 462a + 512a 2 + 144a 3 )C 2 + 144(6a - 1)C 3 } 
+12Ci{(81 - 190a + 200a 2 + 80a 3 )C 2 + 24(7a - 1)C 3 } , 

= 9(61 - 162a + 100a 2 + 52a 3 )(j (3.7) 
+4C 2 {(185 - 414a + 352a 2 + 144a 3 )C 2 + 288(3a - 1)C 3 } 
+12Ci{(125 - 322a + 256a 2 + 64a 3 )C 2 + 12(22a - 7)C 3 } , 

= 3(7 - 16a + 16a 2 - 4a 3 )C 2 + 4(2a - 1)C| + 4CiC 2 (2a + 1) . 
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Unfortunately, the only real set of (Ci>C2>C3>«) that solves (3.7) is 

C 2 = -1.03Ci, C3 = 2.01Ci, a = 2.65, (3.8) 

and we get m 2 < 0, i.e., a tachyon instability. Therefore nPF theory with a generic polynomial 
3-graviton interaction is inconsistent either due to nonrenormalizability or instability. 

3.2 Derivative + polynomial interaction 

A more general interaction may contain derivatives, and a reasonable way to introduce such 
interactions is to expand ^/—gR to higher orders in h with g = i] + jj^. That is, we now 
investigate 

S = S nPF + 2M 2 J d 4 x[[^R]\ h3 + [^R]\ h4 

2 
2 

-^{yiKKKK + ite(V^) 2 + vMKKh + y A h^h 2 + y ^)] 

= SnPF + J d 4 x (iV 3h hhh + iV 4h hhhh) . (3.9) 

Now that Mp l plays the role of the coupling, we need the cubic and quartic vertices in order 
to get loops of 0(Mp 2 ). Explicit form of y/—g R\, 3 and \J—g i?L 4 can be found in, e.g., [15], 
and and V 4 /j as well as ingredients of the 1PI 1-point and 2-point functions are given in 
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A. 2. Then the condition (3.2) is 

= 536 + 333x 2 + 36xi(54 + 27x 2 - 8x 3 ) + 4x 2 (402 + 193x 2 - 144x 3 ) + 1032x 3 

-a{2104 + 999x1 + 4x 2 (1664 + 655x 2 ) 

+12xi (638 + 271x 2 - 192x 3 ) + 48(97 - 84x 2 )x 3 } 

+2a 2 {l586 + 5694xi + 567x 2 + 5092x 2 + 2340xix 2 
+1948X 2 . + 48(58 - 21x x - 36x 2 )x 3 } 

-4a 3 {559 - 72x1 + 30xi(65 + 12x 2 ) + 4x 2 (439 + 92x 2 )} 

+a 4 {632 - 756x1 + 48xi(51 - 20x 2 ) + 64x 2 (43 - 9x 2 )} , 
= 138 + 549x 2 + 6a?i (145 + 250x 2 - 168x 3 ) + 740x^ + 16x 2 (67 - 72x 3 ) + 156x 3 

-a{598 + 2007xi + 2396x^ + 12x x + 960x 3 + 48x 2 (95 - 96x 3 )} (3.10) 

+2a 2 {556 + 1179x 2 + 6xi(433 + 578x 2 - 264x 3 ) 
+4x 2 (832 + 383x 2 - 432x 3 ) + 888x 3 } 

-4a 3 {253 + 108x1 + 4x 2 (233 + 52x 2 ) + 6xi(131 + 96x 2 )} 

+4a 4 {90 + 248x 2 + 3xi(68 - 39x x ) - 48x 2 (4x x + 3x 2 )} , 
= 28 + 63x 2 + 4x 2 (2 - 3x 2 ) + 12x x (7 + x 2 ) 

-o{12(2 + 3xi)(3 + 4xi) - 8x 2 (5 + 3xi) - 24x^} 

+12a 2 (3 + 14xi + 12x 2 - 4x 2 ) + a 3 (8 - 48x x - 36x 2 ) , 

which admits numerous solutions, such as 

xi = 1.00, x 2 = -4.06, x 3 = 1.68, a = 0.85 . (3.11) 
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For this choice of x's and a, the tadpole condition (3.6), 

= 65 + 3xi (191 + 204xi + 96x 2 - 2808x 3 ) + x 2 (822 - 816x 2 - 12528x 3 ) 

-108x3(23 + 36x 3 ) + 216yi - 144y 2 + 1512y 3 + 2304y 4 + 1728y 5 
-a{472 + 3xi (1483 + 2442x x + 4968x 2 - 10368x 3 ) 

+2x 2 (5145 + 2988x 2 - 25920x 3 ) - 108x 3 (73 + 288x 3 ) 

-408yi - 6000y 2 + 12096?/ 3 + 18720y 4 + 15552y 5 } 
+a 2 {l273 + 12xi (751 + 2010xi + 4974x 2 - 2916x 3 ) 

+12x 2 (3029 + 2948x 2 - 5616x 3 ) - 648x 3 (ll + 96x 3 ) 

-12624yi - 43104y 2 + 33048y 3 + 50976y 4 + 41472y 5 } 
-2a 3 {l283 + 6xi(355 + 2946xi + 7512x 2 - 864x 3 ) (3.12) 

+24x 2 (1243 + 1258x 2 - 288x 3 ) - 1728x 3 

-21552yi - 59520y 2 + 21168y 3 + 30528?/ 4 + 13824y 5 } 
+4a 4 {llll - 3xi(443 - 2307xi - 4896x 2 + 1296x 3 ) 

+12x 2 (1093 + 780x 2 - 432x 3 ) - 432x 3 

-15600yi - 39840y 2 + 7560y 3 + 10080y 4 } 
-8a 5 {539 - xi(903 - 1908xi - 3168x 2 ) + 12x 2 (241 + 160x 2 ) 

-5736yi - 13920y 2 + 1296y 3 + 1728y 4 } 
+32a 6 {49 - 3xi (29 - 12x x ) + 108x 2 - 456yi - 1056y 2 } , 

becomes 

= -3407 + 197yi + 398y 2 + 293y 3 + 639y 4 + 1492y 5 , (3.13) 

and it can be solved for various choices of y's. Since the fa — > ff decay channel is still available 
in (3.9) and a from (3.11) lies within (|, l) , the LW mechanism is applicable. Therefore, nPF 
with both derivative and polynomial interactions may give a unitary and consistent theory 
of an effective quantum massive gravity. 

4 Discussion 

The solidity of the LW mechanism has not been fully established yet. Although there is no 
known example against the LW mechanism, and the consistency of several specific models, 
e.g., [13], was thoroughly verified, other cases, including nPF, need to be worked out in this 
respect. 

Once its consistency is checked, nPF might deserve more attention, because it has ad- 
vantages over other dark energy or modified gravity models. 

• Short distance behavior: To modify gravity at large distances, an extra scalar DOF 
is a requisite. But once incorporated into the gravity sector, the universal coupling 
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of this scalar to sources persists at all scales, ruining the matching of the theory to 
the observations and experiments at the solar system scale and below. Therefore, such 
models need to invoke external mechanisms such as Vainshtein[16] or chameleon [17], to 
make themselves viable. But as mentioned in the introduction, we can obtain a smooth 
limit to GR automatically in nPF: Although nPF has two extra scalars, one normal 
and the other ghost, at distances much smaller than m" 1 the effects coming from the 
coupling of the ghost to matter cancel those from the normal scalar-to-matter coupling. 

• Radiative stability: Generically any dark energy or modified gravity model requires a 
small scale to be introduced to explain the cosmic acceleration at the Hubble scale. 
Keeping such a small scale safe from radiative corrections from, e.g., the Standard 
Model fields is nontrivial and sometimes needs a fine-tuning at the same level of the 
cosmological constant problem. On the other hand, the loop corrections [4] to any gravity 
theory from the Standard Model fields minimally coupled to the graviton take the form 

of r 

~ / d 4 x v /3 5(A + ai R + a 2 R 2 + a 3 R^R» u + •••)• (4.1) 

Then the renormalized cosmological constant can get large due to A in (4.1), but in 
a massive gravity theory the Yukawa suppressed gravitational force dilutes it at large 
distances, solving the cosmological constant problem. At the same time, (4.1) does not 
contribute to the nPF graviton mass term 

J d 4 *v^( - ^) (V^ " ah 2 ) , (4.2) 

for a / |, because the contribution to h iiv h yw or h 2 from the perturbative expansion of 
(4.1) is only 

J dW=M(v^-^ 2 )- ( 4 - 3 ) 

In order for the LW mechanism to work, we require | < a < 1, so that m does not get 
affected by (4.1). In other words, since the matter loop corrections should be covariant, 
they cannot renormalize the non-covariant entity, such as the graviton mass, so that 
m is radiatively stable under the minimal gravity-to-matter coupling. But we cannot 
claim complete radiative stability of m yet: We presented nPF only in a perturbative 
form, and a consistent, nonlinear completion of nPF may require a certain nonminimal 
coupling of gravity to matter, with which m could get radiative corrections. 

With these benefits and peculiarities, nPF might be able to play a role not only in cosmology, 
but in solar system anomalies such as the Pioneer anomaly and even in the Standard Model- 
gravity interactions. 

The analysis so far is performed on a Minkowski background as a proof of concept. But 
in order to apply it to cosmology, we need to generalize it to curved backgrounds. This gen- 
eralization seems straightforward, although the actual algebra gets complicated since we can 
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no longer work in the 4-momentum space; The theory will be defined with the covariantized 
version of (1.1), and we may expect that the decay of the ghost through the LW mechanism 
will be seen by an imaginary shift of the ghost mass due to renormalization. There is also 
an advantage in working in a curved spacetime. For example, if we work in (Anti-)de Sitter 
background, a bare theory has one more parameter, the (A)dS curvature, so that the require- 
ments for renormalizability, (3.2) and (3.6), may be relaxed. But how the LW mechanism 
will work in a curved background needs more investigation, because there are subtleties in 
quantum field theory in curved spacetimes, e.g., the preparation of a vacuum. It would also 
be interesting to look into whether nPF can evade generic problems of Pauli-Fierz massive 
gravity in curved backgrounds, such as the Higuchi ghost [18] or some issues raised in [19]. 
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A nPF massive graviton loops 

We can add arbitrary cubic and quartic interactions to S'dpf by introducing the cubic and 
quartic vertices, V 3 ^ and V^: 



S = S„p F + J d 4 x (i\~ 3h hhh + iV Ah hhhh) . (A.l) 



To find loop corrections to SnPF, we path-integrate over h, with 0(h 3 ) and 0{h A ) terms 
providing quantum interactions: 



Z[J] = Jvhexp i(s + J d^r/i, 

= N exp ^ i J d 4 z{*V 3 , (-^y) 3 + *V 4 , 



x exp 



% - J d 4 xd V(x) ( - iP(x - y)) J(y)] , (A.2 



L2 

where W is a normalization constant and P is the inverse Fourier transform of the tree level 
nPF graviton propagator, (1.2). With 

P = x/x>w, ij d 4 xJ(x)( - iP(x - y)) = Xn^% , (A.3) 

we can expand (A.2) diagrammatically. The pieces of our interest are 



-6 (J --(l8 +18xx~<[ V— x|, (A.4) 



(a) (b) (c) (d) 

where (a) ~ (d) are the loop parts of corresponding diagrams. 
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A.l Polynomial interaction 

For the interaction of (2.6), there is only the 3-point vertex 

= _ (2!) 3 3! (Cl 7 ?A t 1^3 r /^lA t 2 r /^2A 1 3 + C2 7 ?/ilA12 7 ?l y 1^2 7 ?A t 3^3 + C3 ? ?/.ti^i 7 ?/i2l'2 7 ?A t 3l'3 

+symmetrization in \iv + permutation in 123) , (A. 5) 

and we have only to consider diagrams (a), (c) and (d). Among their divergences, only that 
of (0?) has 0(k 2 ) pieces: 

1 ?'A 2 r 

{27(9 - 30a + 28a 2 + 4a 3 )C 2 



e 69127r 2 m 2 (l - a) 3 

+4C 2 ((179 - 466a + 512a 2 + 144a 3 )C 2 + 144(6a - 1)£ 3 ) 

+12Ci((67 - 186a + 168a 2 + 80a 3 )C 2 + 24(7a - l)( 3 )}k 2 ^ 

+{9(5 + 8a - 16a 2 + 36a 3 )C 2 + 12(7 - 2a + 16a 2 )CiC2 + 4(7 + 2a)C|}A: 2 y 
+{9(15 - 44a + 48a 2 - 28a 3 )C 2 - 12(1 - 6a + 8a 2 )CiC2 - 4(11 - 14a)C|}l 3 
+ {9(17 - 36a + 8a 2 + 20a 3 )Ci 2 + 4 C2((3 + 26a - 80a 2 )C 2 - 72C 3 ) 

+12Ci((29 - 68a + 44a 2 -8a 3 )C 2 + 6(8a-5)C 3 )}l4 , (A.6) 
from which we can read off b\ 6 2 6 3 and 64. 
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A. 2 Derivative + polynomial interaction 

For the more complicated interaction of (3.9), we have 

~ 2^' ll/ 2^2/li ? ?^i/^2 r //J3l / 3 "I" 2^ 1 ^3^2l/ 3 ? //Ji/X2 r ?^l^2 "I" ^l/i3 ^VsVuiHaVflVS (A-7) 

~\~2kl^ 2 ^21/3 ^11/2 ^1/^3 "I" ^1/X2 "bfllVvi V3 r lv2^2, ~ ^1 ' ^2^1/^2 Vu2P3 Vv3fJ,l 
,2 



m 
~~2 



+symmetrization in /itz/ + permutation in 123 , 
i f 1 1 

^4ft = ~~ J^2~\ ~~ lg^ 1 ' ^ 2 ^Ml !y l 7 ?/ i 2^2 r //i3t / 3 r ?M4l / 4 ~~ g^lM2^1l'2 r //ill'l 7 ?M3^3 7 ?M4^4 

~\~~gki ■ k2f]^ ll/1 i]^ 2U2 i]^ 3 f li rj U3Ui + -^k\ lJL2 kiy 2 rj^ lVl rj^ 3l j lA rj V3UA 

Pui^viv-i^puiv/i, + r) kip l2 ki l , 2 r] l j il ^ 3 r] 1 y lU3 T]^ iUi 

~\"^k\^ 3 k2ii 3 Tj^ 1 ^ 2 i]y lU2 i]^ AU4 — —ki ■ k2f] v - L fj, 2 rj V2IX3 rj V3 fj il ri i j iAVA 

1 1 

~\-~ki^ 2 k2p Jl Tji /2 ^ 3 rj L , 3Ul r]^ 4 ^ 4 + — ki^kip-s'n^Ki^piviV^^ (A. 

— fei • k2f]p LXVx T\ V2 p L3 T\ V3 p l ,/j} V4il2 — 2ki )12 k2[ 13 T] U2 n 4 r} U4fll ri UlU3 

~^2^ 1 ' ^^iMS^^^/M 7 /^/^ ~~ O ' ^ 2r7 A t 1^3 r ?^1^3 r ?M2At4 T /^2^4 

~2^ l ^s^P4 r l^a^4 r lpip,2 r luiU2 ~ 2ki fl2 k2 f i 3 'r] U2f j jl ri ulll4 ri UiU3 
-\-k\ • k2 r ri Ul ^ 2 rj U2 ^ 3 rj V3l j Li rj Vi ^ 1 — kxij, 2 k\ V2 T] uli j i3 ri U3lXi r] UAIX1 

— k\y i3 k2v 3 y}v\li/J)v^li2 r \v2p,\ + klfJ,2^P3 r t^2Pi r JuiU 3 r JfJ.4U4 ~ ^l/x 3 ^2/x 4 ^3^2 V^Pl Vvi V4 

—kifj, 2 kif J , 3 Tiij, 1 i, 1 r] U2 ij lA Ti UAl/3 — k\^ 2 k\i 13 ri l , 2l , 3 ri^ ll j Li rj VlVA — ki^ 2 k2^ 1 rj l/2 ^ 3 ri l/3 ^ 4 ri UiUl 

2 

m , 

2~ (,yi r ?MlM2 ? 7l'lM3 7 7! / 2A'4 7 ?i / 3^4 "I" y2Vp,lP2 7 ]viU2 r lp3lJ.4 r lv3V4 + 2/3 ? ?A i l/i2 ^1/^3 ^2 ^3 ^4 ^4 



+2/4 7 ?/ii/^2 7 /l / 1^2 7 /A t 3^3 7 /A t 4^4 ~l~ y^> 1 lpiUl 1 lp,2f2 r lp,3U3 r IP4^4) 

+symmetrization in + permutation in 1234| . 
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Then the divergences up to 0(k 4 ) are 



1 im 



0) = 977 n + 21x i + 30x 2 + 12x 3 - a(23 + 66x 4 + 96x 2 + 48x 3 ) 

e 2ob7r z Mp(l — a)- 3 L 

+4a 2 (7 + 18xi + 24x 2 ) - 12a 3 (l + 3xi + 4x 2 )}^ + finite , (A.9) 



im 



2 



e 5127r 2 AfJ, (1 



(6) = 1 (_7 + 22a - 24a 2 + 12a 3 ) f A; 2 Ii - k 2 \ + \ - I 



2 2 



9 

+— {3(3yi - 2y 2 + 21y 3 + 32y 4 + 24y 5 ) - 2a(14yi - 8y 2 + 99y 3 + 156y 4 + 144y 5 ) 



+4a 2 (2 - a)(4 yi - 4y 2 + 27y 3 + 36y 4 )}li 

— {3(7 - 30 Vl - 68y 2 - 3y 3 - 4y 4 ) - a(69 - 284 yi - 640y 2 - 36y 3 - 48y 4 ) 

+4a 2 (21 - 76yi - 176y 2 ) - 4a 3 (9 - 38y 4 - 88y 2 )}l 2 + 0(k 4 )] + finite , (A.10) 



im 2 



(C) = 1 1536^3(1 -.Wa-l) * 7 + 2111 + 3012 + 1213 " 0(23 + 6611 + + i8X3) 
+4a 2 (7 + 18xi + 24x 2 ) - 12a 3 (l + 3x 4 + 4x 2 )} 

+ finite , 



k% - k 2 ^- + ^ - I 4 ) + 2m 2 (x 2 + 6x 3 )Ii + ^(3xi + 4x 2 )I 2 + 0{k 4 ) 



(A.ll) 
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(d) - l ^mui-a)' K 462 + 243x ? + 1211(33 + 6712 - 2413 

+2(792x2 + 358x1 + 516x 3 - 288x 2 x 3 ) 
-a(l978 + 1053x 2 + 6304x 2 + 2580x1 

+12xi (138 + 253x2 - 192x 3 ) + 4656x 3 - 4032x 2 x 3 ) 
+2a 2 (l592 + 783x 2 + 6x^233 + 354x 2 - 168x 3 ) 

+4(1167x2 + 489x1 + 696x 3 - 432x 2 x 3 )) 
-4a 3 (623 + 162x 2 + 1572x 2 + 368x1 + 6xi(107 + 44x 2 )) 
+4a 4 (206 - 27x 2 + 688x 2 - 144x1 + 60xi(7 - 4x 2 ))}fc 2 Ii 



+{37 + 774xi + 45x 2 + 12x 2 + 84xix 2 + 28x1 
-a (63 - 27x 2 + 176x 2 + 20x1 + 12xi(250 + 9x 2 )) 
-a 2 (6 + 216x 2 - 424x 2 + 8x1 - 24xi(179 + 9x 2 )) 

+4a 3 (32 + 117x 2 - 92x 2 - 6xi(109 + 8x 2 )) - 12a 4 (8 - 32x x + 27x 2 )}/c 2 I 2 
+{103 + 270x 2 + 28x 2 - 88x1 - 6xi(59 + 4x 2 ) 

-o(437 + 1062x 2 - 336x 2 - 200x1 - 24xi(65 + 7x 2 )) 
+2a 2 (273 + 828x 2 - 432x 2 - 56x1 - 42x x (29 + 4x 2 )) 

-4a 3 (67 + 342x 2 - 152x 2 - 48xi(8 + x 2 )) +8a 4 (7- 18xi + 63x 2 )}l 3 (A.12) 
-2{162 - 153x 2 + 256x 2 - 12x1 - 3x x (79 + 116x 2 - 120x 3 ) + 438x 3 + 288x 2 x 3 

-a(690 - 477x 2 - 12xi(80 + 97x 2 - 78x 3 ) + 4(218x 2 + 23x1 + 462x 3 + 72x 2 x 3 )) 

+4a 2 (259 - 99x 2 + 335x 2 + 106x1 - 12xi(25 + 28x 2 - 12x 3 ) + 474x 3 ) 

-4a 3 (185 + 27x 2 - 12x x (6 + 13x 2 ) + 80x 2 (4 + x 2 )) 

+4a 4 (58 + 45x 2 + 220x 2 + 12x x (9 - 2x 2 )) }l 4 
+2{l + 306x 2 + 528x2 + 672x1 + 216x 3 + 648x 2 x 3 + 648X 2 . + 12xi(16 + 75x 2 + 27x 3 ) 

-3a(l0 + 405x 2 + 12xi(24 + 97x 2 + 36x 3 ) + 4(190x 2 + 197x1 + 90x 3 - 432x|)) 

+a 2 (l33 + 1818x 2 + 4008x1 + 864x 3 + 72x 2 (49 + 72x 3 ) + 48xi(31 + lllx 2 + 81x 3 )) 

-12a 3 (15 + 99x 2 + 220x 2 + 148x1 + 12xi(9 + 19x 2 )) 

+2a 4 (38 + 261x 2 + 432x2 + 432x1 + 24xi (10 + 27x 2 ))}m 2 Ii 
+2{187 - 45x 2 + 24x 2 - 156x1 + 12xi(31 - 15x 2 ) 

-6a(l25 - 27x\ - 4x 2 - 68x1 + 6xi(39 - 14x 2 )) 

+a 2 (llll - 234x 2 - 144x 2 - 240x1 + 24x x (83 - 12x 2 )) 

-48a 3 (15 - 9x 2 - 2x 2 - 2x1 + 3xi(9 - 2x 2 )) 

+4a 4 (43 + 84xi - 18x 2 )}m 2 I 2 + 0{k 4 ) + finite . 
From (A. 4), we can extract 

B div = -3(a), n div = 12(6) + 18(c) + 18(d) , (A.13) 
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and the corresponding frj's can be obtained straightforwardly. 
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